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Abstract. We review our recent claim that there is evidence of non-Gaussianity in 
the 4 Year COBE DMR data. We present some new results concerning the effect of 
ly-j ' the galactic cut upon the non-Gaussian signal. These findings imply a localization of 

f^i , the non-Gaussian signal on the Northern galactic hemisphere. 
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I EVIDENCE FOR CMB NON-GAUSSIANITY 

In a recent paper [1] we showed that the 4 Year COBE DMR data exhibits 
evidence of non Gaussianity at a high confidence level. We made use of statistical 
tools described in more detail in [2,3]. Since then our result has been corroborated 
by two other groups [4,5]. In this review we revisit our analysis, and the tests to 
which we have subjected it, including some new results. 

In our analysis we propose, and work with, an estimator for the normalized 
bispectrum, denoted by If. We refer the reader to [1] for its definition. We then 
applied this estimator to the inverse noise variance weighted, average maps of the 
53A, 53B, 90A and 90B COBE-DMK channels, with monopole and dipole removed, 
at resolution 6, in ecliptic pixelization. We use the extended galactic cut of [6], and 
[7] to remove most of the emission from the plane of the Galaxy. We apply our 
statistics to the DMR maps before and after correction for the plausible diffuse 
foreground emission outside the galactic plane as described in [8], and [9]. 

By means of Monte Carlo simulations we also found the distributions Pi(If) for 
what we should have seen assuming a Gaussian signal, which is then processed 
by the experimental set up associated with DMR. These P{If) were inferred from 
25000 realizations (see Fig. 1). The observed If and the distributions Pi(If) are 
plotted in Fig. 1. One immediately notices the presence of a significant deviant. 
In order to quantify this deviant we define the goodness of fit statistic 
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FIGURE 1. The vertical thick dashed line represents the value of the observed if. The solid 
line is the probability distribution function of if for a Gaussian sky with extended galactic cut 
and DMR noise. 

X 2 = = l^(-2 log P e (lf) +/%), (1) 

where the constants Pi are defined so that for each term of the sum (Xf) = 1. The 
definition reduces to the usual chi squared for Gaussian Pg. We build a X 2 for the 
COBE-DMR, data from the Pg(Ig) inferred from Monte Carlo simulations, taking 
special care with the numerical evaluation of the constants Pg. We call this function 
Xcobe- We then find its distribution F{X"^ OBE ) from 10000 random realizations. 
This is very well approximated by a x 2 distribution with 12 degrees of freedom. We 
then compute Xq OBE with the actual observations and find Xq OBE = 1.81. One 
can compute P{Xq OBE < 1.81) = 0.98. Hence, it would appear that we can reject 
Gaussianity at the 98% confidence level. 

II IS IT A SYSTEMATIC EFFECT? 

We checked that this result could not be due to the following systematics: 

1. Foregrounds contamination: 

• Dust (using the DIRBE sky maps and also the Schlegel et al dust 
model) 

• Synchrotron (with the Haslam template) 

• Foreground corrected maps (effect persists on corrected maps) 

2. Noise model: 

• Anisotropic sky coverage 



• Noise correlations between different pixels 

• Analysis of noise templates 

3. Galactic cut: 

• Dependence on shape ( "custom" versus constant elevation) 

• Dependence on elevation 

• Dependence on monopole and dipole subtraction, before or after 
the cut, with or with out galaxy. 

4. Possible small residual errors in corrections for 

• Spurious offsets induced by the cut. 

• Instrument susceptibility to the Earth magnetic field. 

• Callibration errors . 

• Errors due to incorrect removal of the COBE Doppler and Earth 
Doppler signals. 

• Errors in correcting for emissions from the Earth, and eclipse 
effects. 

• Artifacts due to uncertainty in the correction for the correlation 
created by the low-pass filter on the lock-in amplifiers (LIA) on 
each radiometer 

• Errors due to emissions from the moon, and the planets. 

5. Assumptions in Monte Carlos: 

• Dependence on power spectrum tilt 

• Dependence on smooth versus discontinuous power spectrum 

• Dependence on beam shape 

• Dependence on pixelization. 

In fact, the confidence level quoted above reflects the worse line up of systematics. 
If we try to correct for systematics, in general the confidence level for rejecting 
Gaussianity is enlarged to beyond 99%, as we describe in more detail in [3] 

III WHERE IS THE NON-GAUSSIANITY? 

We now concentrate on a subset of tests involving the galactic cut which we 
applied to our result. Changing the galactic cut affects sample variance besides 
eliminating possible contaminations from the map. We considered extensive varia- 
tions of the cut, including additions of polar cap cuts to the extended cut. 

We found that cuts from the pole affect the result more than cuts from the 
equator. This suggests that the effect may be localized near the Poles. We therefore 
decided to compare the effect of applying cuts only in the North or South galactic 
poles. We considered cuts down to 60° (2668 pixels excluded). We find the curious 
result that cutting Northern caps is more damaging for the nonGaussian spike than 
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FIGURE 2. The effect of combining the extended cut with polar cap cuts down to 60° in the 
North or South hemispheres. 

cutting Southern caps (fig. 2). Indeed the first few Southern cap cuts appear to 
increase the spike. 

Non-Gaussianity could in principle be localized in Fourier space without being 
localized in real space 1 . Examples of such behaviour are given in [10]. We believe 
that the signal we have found is essentially localized in Fourier space. However the 
results we have just presented suggest that our signal may indeed be also localised 
in real space, around the Northern galactic cap. 
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